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O ! Abstract 

■ The 7] — ► 7r°77 decay was investigated by the SND detector at VEPP-2M e + e~ 



collider in the reaction e + e~ 777. Here we present the results and some details 

of this study. We report an upper limit (90% c.l.) Br(r/ — ► 7r°77) < 8.4 x 10~ 4 as our 
*^ final result. Our upper limit does not contradict the earlier measurement by GAMS 

spectrometer. To facilitate future studies a rather detailed review of the problem is 
also given. 

1 Introduction 

The experimental history of the r\ — > 7r°77 decay begins in 1966 when 
1.2 GeV/c momentum tt~ beam from CERN proton synchrotron was used 
to produce r] mesons on a hydrogen target via charge exchange reaction 



7r + p — > rj + n. (1) 

rj meson was tagged by means of neutron time-of- flight measurement. The 
7-ray emission spectrum from the subsequent r\ decays was measured by a 
lead glass Cherenkov counter. As it was claimed in Ref.[|Tj, the observed 
7-ray spectrum could not be explained by j] — > 27 and 77 — ► 3ir° decay 
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modes only. An equally intensive 77 — ► 7r°77 decay was required to get a 
good approximation of the measured spectrum. 

The next experiment carried out in Brookhaven [j2j utilized a similar 
technique except for an order of magnitude higher incident ir~ beam mo- 
mentum and a spark chamber as a 7-ray detector. No 77 — > 7r°77 signal was 
observed in the collected 47 events. 

That was the beginning of the 77 — > 7r°77 decay puzzle which lasted 
for the next 15 years. Several other experiments were performed during 
this time [3 -r- 13], but the situation still remained unclear. The results 
of these experiments are shown in Fig.|T|. For experiments that did not 
measure this branching ratio directly we calculated it using their results 
and current PDG table values for Br(r] -> 2 7 ) = (39.21 ± 0.34)% and 
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Figure 1: The history of the 77 — > 7r°77 decay measurements up to 1980. 
In many experiments the reaction ([]]) with rather low 7r~ momentum 



~ 1 GeV/c was used as a source of 77 mesons and optical spark chambers 
as 7-ray detectors, as in |2j]. Rather low detection efficiency for photons 
lead to a severe background from intense 77 — > 3ir° decay. This was clearly 
demonstrated in [|TTj. It is interesting to note, however, that despite the 
authors' firm conclusion that all observed 47 events in the 77 peak were 
likely from the 77 — > 3tt° decay they presented an odd r(?7 — ► ^77) = (1.6 + 
4.7)%r(?7 — ► neutrals) as their final result rather than just an upper limit. 
We mention this curious fact here only to show how great is sometimes the 
pressure of previously published results — a source of biases in high energy 
physics experiments not always appreciated. 

High 7 detection efficiency provided by heavy liquid bubble chamber 
techniques was also used in 77 — ► 7r°77 decay studies. For example, for a 
xenon bubble chamber an average 7 detection efficiency was as high as 94% 
T2f1 - In this experiment, the chamber was irradiated by a 2.34 GeV/c ir + 



beam and 77 mesons were produced on quasi-free neutrons in the xenon 
nuclei through the reaction tt + + n — ► 77+73. Experiment |T3j also used 
a xenon bubble chamber but tt~ + 77,^ tt~ + 72 + 77 reaction to produce 
77-s by an incident 3.5 GeV/c Tt~ beam. A variation of this technique, a 
hydrogen target inside a heavy freon (CFsBr) bubble chamber, was used in 
|5|, |7j with the 77 meson production reaction tt + + p — > ir + + p + 77. It was 
again emphasized in [J7J that the experimental difficulty in the 77 — > 7r°77 
decay studies comes essentially from the 77 — > 3tt° background. 

Let us mention also two other bubble chamber experiments, for com- 
pleteness. In H, 0.7+0.8 GeV/c K beam at BNL Alternate Gradient Syn- 
chrotron (AGS) facility was used to produce 77 mesons in a hydrogen bubble 
chamber via reaction K~ +p — > A + 77. In [|J the 7r + + d — > 73+73 + 77 reaction 
in a bubble chamber filled with liquid deuterium was used, 0.82 GeV/c tt + 
beam provided also by AGS. 

These bubble chamber experiments also produced controversial results. 
A real breakthrough in the field occurred in 1981 with GAMS experiment at 
Serpukhov [|I4j1 . The 77 mesons were produced again via the charge exchange 
reaction ([[]). But now the incident Tt~ beam momentum was 38 GeV/c. 
This incident energy increase dramatically improved the 77 — ► 3tt° back- 
ground suppression. In previous experiments low energy (~ 1 GeV) 77-s 
were used, which means relatively high thresholds for 7-ray detection. So 
two softest 7-rays could easily escape unnoticed, and this is one possible 
mechanism for 77 — ► 37r° decay to mimic the desired 77 — ► 7r°77 decay. With 
increasing 77 energies, relative threshold for 7-ray detection gets lower, and 
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so does the ij — ► 37r° background. There is also another benefit of higher 
energies: higher invariant mass resolution of the detector for higher photon 
energies significantly improves suppression of the background. An order 
of magnitude lower upper limit Br(?7 — > 7r°77) < 0.3% (at 90% confidence 
level) was reported by this experiment. 

In the next run of the same GAMS experiment a lower 7r~ momentum 
30 GeV/c was used. But because of improvements in the detector and 
accelerator performance significantly higher statistics was collected: about 
6 x 10 5 7] mesons produced via the charge exchange reaction ([I]). At last the 
r\ — > 7r°77 decay was discovered with certainty QT5j . In the mass spectrum 



of the tt°77 system a narrow r] peak was observed over a smooth background 
caused by the 77 — ► 37r° decay. The statistical significance of the effect was 
greater than seven standard deviations. The measured branching ratio was 

Brfa -> 7T°77) = (9.5 ± 2.3) x 10~ 4 . 

After about two years the same experimental statistics was reanalyzed with 
improved reconstruction program for 7-rays from the showers in the GAMS 
spectrometer (jT6j . The improvement provided better treatment of overlap- 



ping showers and allowed to reduce the background further. The new result 
was gg, |T7j 

Br(7/ -> tt°77) = (7.1 ± 1.4) x 10~ 4 . (2) 

Now some words on history of theoretical studies of the 77 — > 7r°77 decay. 
To our knowledge, the first estimate of this decay rate was given by Okubo 
and Sakita in [jTSj . They assumed that the following effective coupling 



4// = «/V (3) 

was responsible for the decay. To estimate the value of the coupling con- 
stant £, they further assumed that the interaction (|3|) was also the source of 
the 77-71"° transition mass through the tadpole diagram which one obtains by 
contracting photon lines in the 777r°77 interaction vertex. The 77 — 71"° transi- 
tion mass itself was estimated by using the unitary symmetry inspired mass 
formula. This approach was hampered by the fact that the tadpole dia- 
gram is badly divergent and one needs some cut-off to calculate its strongly 
cut-off dependent value. By taking the cut-off at the nucleon mass, they 
obtained T(r] — ► 7r°77) ~ 8 eV. But the authors admitted that this esti- 
mate may easily be wrong by an order of magnitude. For comparison, the 
experimental result (0) implies T(r] — ► ^77) = 0.84 ± 0.18 eV. 
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After some years, it was shown in |T9j that vector meson dominance 
model (VDM) with trilinear meson interactions of the VVP and VPP type, 
previously successfully applied to the 77 — ► 7r + 7r~7 and 77 — > 77 decays [|2D| 



predicts too small (~ 10~ 3 ) ratio of the 77 — > tt°77 and 77 — > 77 decay 
widths, indicating r(r/ — > 7r°77) ~ 0.5 eV. So it became clear that a large 
rate for the 77 — ► 7r°77 decay, which was reported by several experiments at 
that time, caused some theoretical embarrassment. As a possible way out 
Singer suggested in |2]J a rather strong quadrilinear WPP meson interac- 
tion which can be dominating in the 77 — > 7r°77 decay. The VDM model was 
further examined by Oppo and Oneda in J22J. It was confirmed that this 
model predicts a small 77 — ► 7r°77 decay width in the 0.3-0.6 eV range. But 
the authors also indicated a possible additional contribution coming from a 
scalar meson intermediate state. There was very little experimental knowl- 
edge about scalar mesons at that time (even nowadays scalar mesons still 
remain the most experimentally unsettled question in the low energy phe- 
nomenology). Using chiral symmetry arguments Gounaris |23j estimated 
scalar meson contribution and obtained r('77 — ^ 7r°77) = 1.0 ± 0.2 eV in 
excellent agreement with the present experimental result (@). But this cor- 
rect experimental result then had years ahead to come. So the Gounaris' 
result was considered as a failure of the theory rather than success. Some 
theoretical work followed |Z|, |25|, £5j which seemed successful in producing 
large enough decay width by using essentially scalar meson contribution. 
The new experimental situation [[T3|, |I5fl , however, promptly depreciated all 
such attempts and stopped them forever. 

Meanwhile it became increasingly evident that the chiral symmetry and 
its spontaneous breaking play crucial role in the low energy hadron phe- 
nomenology. This was culminated by creating a theoretical framework 
called chiral perturbation theory [^7f where the degrees of freedom used 
are the low-lying hadron states instead of useless in the low energy world 
(current) quarks and gluons. The situation here resembles [|Z8f theory of 
superconductivity: nobody questions the QED as underlying theory of all 
electromagnetic phenomena, but it does not help much while dealing with 
superconductivity. The reason is that electrons and photons are not ade- 
quate degrees of freedom for the superconductive phase. The translational 
symmetry is spontaneously broken due to presence of atomic lattice. The 
resulting Goldstone bosons (phonons) mediate the most important interac- 
tions. Electron gets a dynamical mass and becomes quite different from the 
electron of the bare QED Lagrangian. In the case of chiral perturbation 
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theory (ChPT), the Goldstone bosons, associated with the spontaneous 
breaking of the chiral symmetry, are pseudoscalar mesons. So their inter- 
actions can be predicted to a great deal in the low energy limit. Therefore 
it is not surprising that the last chapter of our theoretical history is related 
to ChPT. 

But at first let us mention some before-ChPT-age papers, for complete- 
ness. It was shown in [j29|, f5tj that, unlikely to 77 — > 77 decay, baryon loop 
contribution vanishes for the 77 — > 7r°77 decay in the chiral symmetry limit. 
Pion loop contribution turned out to be very small |3{J . Failing to find any 
large contribution in the 77 — > 7t°77 decay rate in otherwise successful chiral 
model, the authors of suggested that something may be wrong with the 
experiment itself and soon the GAMS experiment showed that they were 
correct. Right after this crucial experiment two papers appeared [j3T], EEZfl 
which calculated the 77 — > 7r°77 decay width by using (constituent) quark 
loops to estimate various meson amplitudes. They obtained the correct 
magnitude of the decay width ~ 1 eV. 

In ChPT the 77 — > 7r°77 decay was studied in |33j . This decay is rather 
peculiar from the point of view of the chiral perturbation theory. In the 
momentum expansion, the leading 0(p 2 ) term is absent because there is no 
direct coupling of photons to 7r° and 77. Due to the same reason the 0(p A ) 
tree contribution is also absent. So for the 77 — > 7r°77 decay ChPT series 
starts with 0(p 4 ) one- loop contributions. But these loop contributions are 
also very small. Usually dominant pion loops are suppressed here because 
they contain G-parity violating r]7r + 7r~7r° vertex. Kaon loops, although not 
violating G-parity, also give negligible contribution because of large kaon 
mass. It turned out that the main contribution comes from 0(p 6 ) coun- 
terterms which are needed in ChPT to cancel various divergences. The 
coefficients of these counterterms are not determined by the theory itself 
and should be fixed either from experimental information or by assuming 
that they are saturated by meson resonances (vector meson exchange giving 
the dominant contribution). In |33f the latter approach was adopted which 
gave T(r] ^ 7r°77) ~ 0.18 eV - - too small compared to the experimental 
value (@). But, as was argued in fl33j, keeping momentum dependence in 
the vector meson propagators gives an "all-order" estimate (that is 0(p 6 ) 
and higher) of the corresponding contribution of about 0.31 eV, in agree- 
ment with the old VDM prediction |T9|, ffiZfl . Note that although the kaon 
loop amplitude alone is totally insignificant, its interference with the "all- 
order" VDM amplitude gives still small but noticeable contribution [fM 



6 



There are also scalar and tensor meson contributions (which signs cannot 
be unambiguously determined) and one-loop contribution at 0(p 8 ) which is 
of the same order of magnitude as the 0(p 4 ) contribution, because it does 
not violate G-parity. Taking into account all these mechanisms, the final 
estimate of |33j was 

T(r] 7r°77) = 0.42 ± 0.20 eV. (4) 



Some more details were given in [|35f but the result essentially remained 



unchanged. Contributions of the C-odd axial vector resonances were con- 
sidered in These contributions increase (|j) by about 10%. But this 



is still not enough to reconcile the ChPT result (K) with the experiment. 
Therefore it was suggested in [[77f that the resonance saturation assumption 



may not be valid for 0(p 6 ) chiral Lagrangian (although it works quite well 
at 0(p 4 )). Some indications about validity of this suggestion is provided 
by Bellucci and Bruno in ||ggj| . They calculated the coefficients of the 0(p 6 ) 



counterterms in the framework of the Extended Nambu Jona-Lasinio model 
(ENJL) and obtained 

r(r/^7r°77) = 0.58 ±0.3 eV. 

It is higher than (g) and is in agreement with experiment within one stan- 
dard deviation. However, the same coefficients when calculated by using 
another chiral effective Lagrangian, obtained by a bosonization of the NJL 
model, led |3Sj] to small decay width of the order of 0.1 eV, while the 
77 — > 7T°7r experimental data were reproduced well. Another analysis [f|0 



within the ENJL model framework also produced rather small 0{p Q ) rate 
in the leading 1/N C approximation 

rfa - At) = (0.27JU?) eV. 

On the contrary, Nemoto et al. argued |4T[] that the ENJL model might 
be quite successful in explaining experimental 77 — ► 7r°77 width if one in- 
cludes instanton induced 6-quark interactions, which lead to a strong ?7a(1) 
breaking and nonstandard 77-77' mixing angle close to zero. 

As we have seen, 77 — > 7r°77 decay has dramatic, both experimental and 
theoretical, history. And the story is by no means over yet. From the 
theoretical point of view, this decay provides unique possibility to probe 
higher order effects in ChPT. But at present these effects can not be cal- 
culated without model ambiguities. And the situation will hardly change 
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until more precise and detailed experimental data appear. Such data about 
77 — > 7r°77 decay are also highly desirable for understanding other interest- 
ing rare decays, such as Kl — ► 7r°e + e~ j42j (in which, it is believed, a signal 
of the direct CP-violation can be observed) and 77 — > tt°1 + 1~ |E3|, fH, [|5j 
(which is a sensitive probe of C-conservation in the electromagnetic in- 
teraction of hadrons). Despite considerable experimental effort, only one 
reliable measurement of the 77 — > tt°77 decay rate exists at present. There- 
fore p6| , "a re-measurement of the decay rate and a measurement of the 



decay distributions is certainly desirable" 



2 Detector and experiment 

The SND detector (Fig.g) is a general purpose nonmagnetic detector which 
successfully operates at the VEPP-2M collider [§7| since 1995. The de- 
tailed description of this detector and some of its subsystems can be found 
elsewhere [|9|, |5Uj, and will be not repeated here. But maybe it is 



worthwhile to mention some features of the SND calorimeter, because this 
subsystem is crucial in the analysis like one presented in this article. The 
SND calorimeter consists of 1630 Nal(Tl) crystals with a total weight of 
3.6 tons. Crystals are assembled in three spherical shells with solid angle 
coverage ~ 90% of 47r. Each calorimeter layer includes 520-560 crystals 
of eight different shapes. The total thickness of the calorimeter is 13.5Xq 
(35 cm) of Nal(Tl). 

The calorimeter was calibrated by using cosmic muons [|5Tf and e + e~ — ► 
e + e~ Bhabha events [j52j. The dependence of the calorimeter energy reso- 



lution on photon energy is shown in Fig.|3|. This dependence can be fitted 
as: 

4.2% 



E v ' jE(GeV) 

Several effects influence the energy resolution, such as passive material 
before and inside the calorimeter, leakage of shower energy through the 
calorimeter, electronics instability and calibration inaccuracy, light collec- 
tion nonuniformity over the crystal volume. Reasonably good agreement 
between simulated and measured energy resolutions was achieved after all 
these effects were duly taken into account in the simulation. An average 
energy deposition for photons in the calorimeter is about 93% of their initial 
energy. 
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Figure 2: SND detector, section along the beams: (1) beam pipe, (2) drift chambers, 
(3) scintillation counter, (4) light guides, (5) PMTs, (6) Nal(Tl) crystals, (7) vacuum 
phototriodes, (8) iron absorber, (9) streamer tubes, (10) 1 cm iron plates, (11) scintillation 
counters, (12) and (13) elements of collider magnetic system. 



High granularity of the calorimeter ensures rather good angular resolu- 
tion for photons. A novel method for reconstruction of photon angles was 
suggested in fl55fl - It is based on an empirical finding that the distribution 
function of energy deposition outside the cone with the angle 6 around the 
shower direction has the following form for electromagnetic showers in the 
SND calorimeter (according to simulation) 



E(6) = a-exp(-^6/P) 1 (5) 

where the a and (3 parameters turned out to be practically independent 
of the photon energy over the interval 50 700 MeV. Fig.f| shows the 
calorimeter angular resolution as a function of photon energy. This energy 
dependence can be approximated as 

<r = °; 82 0.63°. 
! !E(GeV) 



High granularity also enables to discriminate merged photons in the 
calorimeter by examining energy deposition profile of the corresponding 
electromagnetic shower in transverse direction. Merging of photons is the 
main mechanism for rj —> 37r° decays to imitate the desired rj — ► 7r°77 
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Figure 3: Calorimeter energy resolution versus photon energy. Energy resolution was 
measured by using e + e~ — > 77 (dots) and e + e~ — ► e + e"7 (circles) reactions. 

process (another one is a loss of photons through openings in the calorime- 
ter). So, for background suppression, it is very important to recognize 
merged photons. In j53j a parameter was suggested to separate hadronic 
and electromagnetic showers in the SND calorimeter by comparison of ac- 
tual transverse energy profile of the shower with the expected one (|5|) for 
a single photon. The same parameter also helps to discriminate merged 
photons. 

The results presented in this article are based on the statistics collected 
by SND in two experiments. The first one was carried out in the time pe- 
riod from February 1996 until January 1997. Seven successive scans were 
performed in the center of mass energy range from 980 to 1044 MeV. Data 
were recorded at 14 different beam energy points. Collected statistics was 
almost tripled during the next experiment in 1997-1998. It consisted of 
three scans at 16 energy points in the interval from 984 to 1060 MeV. The 
total integrated luminosity used in this work is equal to 12.2 pb _1 and cor- 
responds to about 2 x 10 7 (f) mesons produced. The corresponding number 
of 77 mesons from — > 777 decay is 2.58 x 10 5 . The collider luminosity was 
measured using the process e + e _ — ► 27. Systematic error of the luminosity 
measurement estimated using e + e~ — ► e + e _ reaction is close to 3%. 
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Figure 4: Calorimeter angular resolution versus photon energy. 

3 Theoretical models for Monte Carlo simulation 

The general form of the rj(p) — ► 7T°(q)^(ki, 61)7(^2, £2) transition ampli- 
tude as a function of particle momenta and photon polarization vectors is 
dictated by gauge invariance, Bose symmetry, and CP-conservation 



M = ei^e 2v 



A(x 1 ,x 2 )T(^ + 



B{xi,x 2 ) mliV 



-T, 



mi 



(2) 



where A(x\, x 2 ) and B(x\, x 2 ) are Lorentz-invariant form factors, which are 
symmetric functions of dimensionless variables 



X\ = 



p ■ k 



1 



mi 



x 2 



p ■ k 2 



mi 



"1] Tj 

Two independent gauge-invariant tensors are defined as 

Tfi = Kk\ - h ■ k 2 g^, 

Tj£ = h-p klp v + k 2 -p p^k\ -h-k 2 p»p v -h-pk 2 -p . 

If CP is not conserved two more gauge-invariant structures containing 
t\iv\o (pseudo)tensor appear |)55f. However we can safely neglect tiny effects 
associated with additional form factors. 
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The decay width is given by the standard expression 



^ r _ 5(p - q - h - k 2 ) |^| 2 dq dki dk 2 



16 (2tt) 5 m 



E n Ei E 2 



(6) 



It is straightforward to perform summation over photon polarizations and 
obtain 



1 



E\M\ 2 = -m* 

eie 2 - 



A(x h x 2 ) + -B(x h x 2 ] 



m: 



2(xi + x 2 ) + ^-l 



mi 



+ 



1 



B(x h x 2 )\'' 



4xix 2 — 2(xi + x 2 



m 



mi 



£ + 1 



Then the decay width is given by the integral 



rfa - At) = ^ ^ / / ^2 E l-^l 2 , 



where the integration limits are 



Xi — 0. — ~~ I 1 

i ' 1 2 



mf / 2 



i 1 1 — 2a?i 



and 



2(1 - 2xi N 



1 - 2xi 



We need some model to calculate A(xi,x 2 ) and B{x\ 1 x 2 ) form factors 
and proceed further. As was already mentioned in the introduction, vector 
meson exchange r\ — ► V*j — > 7r°77 gives the most significant contribution. 
The corresponding form factors have the following form 



- E 9Vra9v 



1 — X\ 



+ 



1—^2 



71 1 1 - 2*! - 4 1 - 2x 2 - Ml ' ■ 



2~2 9vrn9v 
v=p,u),4> 



Bvdm(%i, x 2 ) — 
1 

+ 



^ 7, l-2x 1 -4 l-2x 2 -4 
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In the model considered in j44j, which we choose as the base for our Monte 
Carlo simulation, scalar meson exchange 77 — ► 7r°ao — > 7r°77 contribution 
was also taken into account. If clq is considered as a point-like source of 
photon emission then only A(x\, Xq) form factor acquires additional contri- 
bution from this mechanism 

A ao (x h x 2 ) = P 9a ^ 



ml — rr& 



2(x 1 + x 2 ) + ^-l 



For the coupling constants we take their central values from Ref. , where 
they were extracted from the experimental information on various decay 
widths. The relative sign of the VDM and ao exchange amplitudes cannot 
be determined this way. We have chosen constructive interference between 
them because the experimental width is larger than VDM prediction. 
The following resulting values were used in the simulation: 

m ^w7SVy = 0.0495, m^g^g*^ = 0.0294, 

m ?5V>7775V7 = 0.00267, g am g ail = 0.0132. 
As a check of the simulation program, we reproduced the result of Ref. []44 
for the decay width T(r] — ► 7V°jj)vDM+a = 0-36 eV. 

Equation (j6|) can be rewritten in a form more convenient for simulation. 
The (5-function (or energy-momentum conservation) enables us to perform 
a trivial integration over dq and d cos 62 (where 62 is the angle between k\ 
and A)2 3-vectors) and get in the ?7-meson rest frame: 

Xi, X2 which were defined earlier, in the ry-meson rest frame coincide with 
photon energies divided by ?7-meson mass. 6\ and (f \ are the first photon 
polar and azimuthal angles in the 77 meson rest frame and (p2 is the sec- 
ond photon azimuthal angle in the plane perpendicular to the first photon 
momentum k\. In fact the integrand in (0) after taking into account the 
energy-momentum conservation depends only on x\ and X2- It is conve- 
nient, however, to introduce another angular variable instead of X2 in the 
following way. Let S* be a rest frame of the system consisting of 7r° and 
the second photon. Of course this system moves in the direction opposite 
to the k\ if initial 77 decays at rest. For the energy E* and invariant mass 
m* of the 7r°7 system we get easily 

E* = 777^(1 — xi), m* = m^yjl — 2x[. (8) 
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Therefore the Lorentz factor of this system is 

7 = = n • (9) 

m* y 1 — 2xi 

The energy of the second photon E\ in the S* reference frame is uniquely 
determined by and hence by X\. 

2m v m* 2y 1 — 2xi 

but its direction is arbitrary. Let us denote the polar and azimuthal angles 
(#2 and (^2) of the second photon in 5* as y>\ = P2 and B\ - - as an angle 
between the velocity 3-vector of the reference frame S* and the direction 
of the second photon. The Lorentz transformation relates X2 and cos Q\ m 
the following way 



X2 = 7^2(1 + /3 cos 6*2), (3= . 

7 

Therefore we can rewrite (0) as 

647T' 3 eie 2 2 27T Z 27T 

where 

dx 



/(xi) = I^¥ 2 = ^ = ^ ~ 1 x ^ 

Having equation QTTf) at hand, we can generate 77 — > 7r°77 decay in the 
7] meson rest frame by the following algorithm: 

• generate first photon normalized energy x\ as a random number uni- 
formly distributed from to xf. Calculate (7r°, 7)-system energy, mass 
and the Lorentz factor according to equations (§) and 



• generate a random number (p* uniformly distributed in the interval 
[0, 2tt] and take it as the azimuthal angle of the S* velocity vector 
in the 77 meson rest frame, which will be called laboratory frame 
below. Generate another uniform random number in the interval 
[—1,1] and take it as a cos/9*, 6* being the polar angle of the S* 
velocity vector in the laboratory frame. This defines the unit vector 
ft = (sin 0* cos tp* , sin 0* sin cp* , cos 0*) along 5* velocity. 
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• now the first photon 4-momentum in the laboratory frame can be 
constructed: E\ = m„xi, k\ = —E\n. 



• generate in the analogous manner ipi, and cos #2, construct the unit 
vector n* along the second photon velocity in the S* frame: 

n* = (sin6 , 2COS(^25 sm ^2 sm( / 9 25 cos ^2)- 

• calculate x\ according to (PDj ) and construct the second photon 4- 
momentum in the S* frame: E% = x^m^, k\ = E^n*. 

• construct the 7r° meson 4-momentum in the 5* frame: E* = m* — 
Eh € = ~kl 

• construct the 7r° meson and second photon 4-momenta in the labora- 
tory frame from their 4-momenta in the frame S* by using appropriate 
Lorentz transformations. 

• for generated x\ and x<i , calculate 

^ = /(^i)EI^I 2 - 

• generate a random number zr uniformly distributed in the interval 
from to z max , where z max is some number majoring z. That is for 
all 77 — ► 7r°77 decay configurations (for all acceptable values of X\ and 
X2) we should have z < z max . In our program we had used z max = 
1.8 x 10" 5 . 

• if z > zr, accept the event, that is the generated 4-momenta of the 7r° 
meson and both photons. Otherwise repeat the whole procedure. 

Note that equation (\n\) indicates that the decay width can be calculated 
by using the above described Monte Carlo algorithm according to 

r(„ - A 7 ) = I ^5 l^l 2 ) (*? - "*")■ 

where < . . . > stands for the mean value. As another check of the program, 
we had reproduced the result of |44j for the width by this method also. 

To understand whether model dependence of the 77 — ► 7r°77 decay matrix 
element can bring a dangerous systematics in our experimental study of this 
decay, we performed another Monte Carlo simulation of this process using 
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significantly different form factors of the quark-box model from [f|5j . These 
form factors are reproduced below (in the units of GeV -2 ): 



2 

Tfl 

2(xi + x 2 ) + —J - 1 



A QB (x h x 2 ) 



-0.616 + 2.14(xi + x 2 ) - 2.509(x? + x\) - 4.184xix 2 + 

1.5896(xf + x\) + 2.936xix 2 (xi + x 2 ), 

Bq B (xu x 2 ) ~ -0.866 + 1.674(xi + ar 2 ) - 3.26(^ + x\)- 

1.781xiX2 + 2.37(x? + x\) + im§xix 2 (xi + x 2 ). 

For the quark-box model, we had used z max = 4.37 x 10~ 5 as majoring 
value in our Monte Carlo program. The prediction of this model for the 
decay width is about two times higher than the VDM prediction and is 
closer to the present experimental result. But this model is also more am- 
biguous. Its predictions strongly depend on the assumed values for the 
u and d constituent quark masses (the above given form factors corre- 
spond to m u = rrid = 300 MeV). It was proven that in the chiral sym- 
metry limit baryon loop contribution strictly vanishes [[29], |30j| . In order 



to check whether quark-box diagaram contribution is indeed dominant we 
must consider other quark-loop diagrams also in the framework of some 
self-consistent model (for example ENJL model) with realistically broken 
chiral symmetry. Therefore we took VDM as a basis for our Monte Carlo 
simulation and used the quark-box model only for systematic error estima- 
tion. 

4 Data analysis and results 

A primary selection of the candidate events for the process 

e + e~ — > (/) — > 777; 77 — ► 7r°77 (12) 
was done according to following criteria: 

• an event must have no charged particles and contain exactly 5 photons 
in the calorimeter; 

• the azimuthal angle of any final photon lies within the interval 27° < 
6 < 153°; 
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• the energy of the next to the most energetic photon is less than 0.8E, 
E being a beam energy; 

• the total energy deposition of final photons is in the range 0.8 < 
E tot /2E < 1.1; 

• the normalized full momentum of the event (Pt t/2E) is less than 0.1; 

The third condition is against QED background e + e~ — ► 27, with extra 
photons originating either from the beam background or splitting of elec- 
tromagnetic showers in the calorimeter. The two latter conditions eliminate 
the main part of the background from the — ► K^Ks decay. 

For the events which passed the primary selection criteria a kinematic fit 
was performed assuming energy-momentum conservation and the presence 
of a 7T° meson in an intermediate state, x 2 °f this fit (x-k) can be used to 
suppress the most dangerous background from 

e + e~ — ► 777; 77 — ► 37T° . (13) 

In particular the Xn < 10 condition being used in our event selection pro- 
cedure reduces the background ( P~3|) by a factor of two. Another ~ 20% 
reduction of this background was achieved by demanding that any photon 
pair not involving photons from the 7r° meson found by the kinematic fit 
has an invariant mass outside m n o ± 30 MeV interval. 

Another significant background left after the selection criteria described 
above comes from the process 

e + e~ -> A; -> 7r°7r°7 . (14) 

This background is only slightly (~ 15%) reduced by the x-n < 10 cut but is 
suppressed by an order of magnitude after the invariant masses of photon 
pairs are restricted as described above to exclude the presence of a second 
7T° meson. 

To further reduce this background, special kinematic fit was applied to 
the selected events assuming utt° intermediate state. Only the events for 
which this kinematic fit fails were selected for further analysis. 

One more background source is the e + e~ — ► 37 reaction (both pure QED 
and 4> — > 777, 77 — > 27 decay) with extra background photons. These extra 
photons are usually soft. To suppress this background, a 37-kinematic 
fit was performed under assumption that two photons in the event were 
spurious. All 37-combinations were examined and one with the best x 2 
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(denoted as X37) was selected. The e + e~ — > 37 background is efficiently 
suppressed by the condition ^37 > 20, while the process ([12]) itself is only 
marginally effected. 

As was already mentioned, transverse energy profile of the shower en- 
ables us to recognize merged photons and additionally suppress the back- 
ground ( |I5|) . The corresponding parameter £ 7 is described in j53j. Fig.|5] 
shows distributions over the parameter £ 7 for simulated signal events flI2|) 
and background flI3p . The efficiency of this parameter in the background 



suppression is clearly demonstrated by this figure. We require £ 7 < — 5 
which reduces background f |I3D by about six times. 




Figure 5: £ 7 distributions. Histogram — simulated background (jT3|) normalized to exper- 
imental statistics. Points with error bars — experiment. Shaded histogram — expected 
rj — ► 7r°77 signal according to simulation of the process flT2|). 



Fig.|5] also shows that almost all background remaining after our cuts 
comes from the process (|I31) . In particular the background from the re- 
action (|P|) seems to become insignificant after application of described 
above selection criteria. But unfortunately it is very difficult to suppress 
the background ( p73|) any further without significant improvement of the 
detector performance (angular and energy resolution for photons). So the 
only option left for us is to subtract the remaining background (fT3| ) relying 
upon the Monte Carlo simulation. 
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For this purpose, special kinematic fit was performed for selected events. 
It was assumed in the fit that four of five photons in each event are from the 
rj — > 7r°77 decay. Of course there are many possibilities how to divide five 
photons into the recoil photon from the — > 777 decay, two photons from 
the subsequent 77 — > 7r°77 decay, and two photons from the 7r° — ► 77 decay. 
The best combination was selected for each event during the fit. x 2 °f the fit 
(Xtt77) is small for desired events from (|I2j) , while background ( |T3j) has wider 
distribution. Fig.0 shows distribution for 170 experimental events left 
after the above described cuts with one more condition Xiry) < 40. The 
simulated background from the process ( |T31 ) has similar distribution over 
this parameter. 



3 
e 

I 12 



18 
16 
UV 



£ 10 




5 10 15 20 25 30 35 40 



Figure 6: x-k^ distribution. Points with error bars — experiment. Shaded histogram 
corresponds to the simulated background (|T3D, normalized to experimental statistics. In- 
sertion shows the distribution for simulated 77 — >■ 7r°77 events. 



Signal and background composition of our data sample was estimated by 
using binned maximum likelihood method incorporating the finiteness of 
the Monte Carlo statistics Using HMCMLL routine from the HBOOK 
package [|57j , the experimental histogram was fitted by the sum of two 
Monte Carlo histograms for the signal (|T2|) and background (|T3j) . As a result 
of this fit, the signal events fraction in the experimental data, depicted in 
Fig.§, was found to be (4.1±7.6)%. This corresponds to the total number of 



19 



signal events of the process ( \T2\ ) N$ = 7.0ig 2 5 9 , where the asymmetric errors 
correspond to the 68.27% confidence interval determined according to Table 
10 from [|58|1 . Number of background events from this fit Nb = 163 ± 19 is 
in a good agreement with the prediction of simulation 158.3 ± 11.5. 

To transform Ns into the branching ratio, we need a detection efficiency 
for signal events. This efficiency determined from simulation of the effect 
with VDM matrix element is e = (14.1 ± 0.4)%. The quark-box matrix 
element gives almost the same value e = (14.8 ± 0.4)%. Therefore the 
systematic error coming from model ambiguity is negligible in comparison 
with the inaccuracy in the background subtraction. The Monte Carlo de- 
tection efficiency was corrected for a loss of events due to additional beam 
background photons causing signal events to be rejected by the selection 
cut requiring exactly 5 photons in an event. This effect not modelled in 
the simulation was as large as 5% in the first SND experiment and 8% in 
the second one with a twice larger statistics. Thus as an average corrected 
detection efficiency for signal events we took e = (13.1 ± 0.8)% with the 
model ambiguity also taken into account. 

The number of 77 mesons in our initial data sample is 2.58 x 10 with a 3% 
systematic error due to limited accuracy of the luminosity determination. 
Given above values for Ns and e we get the branching ratio 

Br(r] -► At) = (2-l±?:|) x 10~ 4 . (15) 

Since the error is very large only an upper limit of the branching ratio can 
be determined from our data: 

Br(n -> 7r°77) < 8.4 x 10~ 4 (90% C.L.), (16) 



where we have again used Table 10 from |5S|, according to which the mea- 
sured mean of 0.54cr implies 90% C.L. interval fi < 2.18cr for the mean \i 
of a Gaussian constrained to be non-negative. 

The systematic error caused mainly by inaccuracy of the simulated £ 7 - 
distribution was estimated to be less than 20% which is negligible in com- 
parison with the statistical error. 

5 Final remarks 

Our investigation once more clearly demonstrates experimental difficulties 
in 77 — ► 7r°77 decay studies caused mainly by the background ([T3|). Our 
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analysis yields the value of the ij — > 7r°77 branching ratio (with systematic 
error not grater than 20%) Br{j] — > 7r 77) = (2. 1 jlf ;|) x 10~ 4 which is con- 
sistent with a more accurate result of GAMS [[16], fT7|l : Br{j] — ► 7r°77) 



7. 1 ± 1.4) x 10 4 , but shows an urgent need of new high precision mea- 
surements. As a final result we present an upper limit: 

Br{j] 7T°77) < 8.4 x 10~ 4 at 90% C.L., 

leaving open the question whether ChPT has problems with this decay. 
Nevertheless, we think that the investigation presented in this paper demon- 
strates feasibility of such studies in future high statistics collider exper- 
iments. Let us remind in this respect that r\ mesons will be copiously 
produced at 0-factories, for example, about 10 8 (f) — > 777 decays per year 
are expected at DA$NE |59[] . We hope that this article will stimulate future 



investigations at ^-factories and will be useful in such studies, 
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